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Abstract 



Operator fractional Brownian motions (OFBMs) are zero mean, operator self-similar 
(o.s.s.), Gaussian processes with stationary increments. They generalize univariate fractional 
Brownian motions to the multivariate context. It is well-known that the so-called symmetry 
group of an o.s.s. process is conjugate to subgroups of the orthogonal group. Moreover, by 
a celebrated result of Hudson and Mason, the set of all exponents of an operator self-similar 
process can be related to the tangent space of its symmetry group. 

In this paper, we revisit and study both the symmetry groups and exponent sets for the 
class of OFBMs based on their spectral domain integral representations. A general description 
of the symmetry groups of OFBMs in terms of subsets of centralizers of the spectral domain 
parameters is provided. OFBMs with symmetry groups of maximal and minimal types are 
studied in any dimension. In particular, it is shown that OFBMs have minimal symmetry 
groups (as thus, unique exponents) in general, in the topological sense. Finer classification 
results of OFBMs, based on the explicit construction of their symmetry groups, are given in 
the lower dimensions 2 and 3. It is also shown that the parametrization of spectral domain 
integral representations are, in a suitable sense, not affected by the multiplicity of exponents, 
whereas the same is not true for time domain integral representations. 

1 Introduction 

This work is about the class of operator fractional Brownian motions (OFBMs). Denoted by 
Bh = {BH{t)}t£M. = {{BnAit), . . . , BH,nit)y G W^,t £ M}, these are multivariate zero mean 
Gaussian processes with stationary increments which are operator self-similar (o.s.s.) with a 
matrix exponent H. Operator self-similarity means that, for any c > 0, 
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{BH{ct)}teR = {c^BH{t)}teR, 




where =c stands for the equality of finite-dimensional distributions and c := e := 
X;^o^''(lnc)'=/A;!. It is also assumed that OFBMs are proper, in the sense that the support 
of the distribution of B}i(t) is M"^ for every t G M. OFBMs play an important role in the analysis 
of multivariate time series, analogous to that of the usual fractional Brownian motion (FBM) in 
the univariate context. They have been studied more systematically by Pitt (1978), Mason and 
Xiao (2002), Bahadaron, Benassi and D§bicki (2003), Lavancier, Philippe and Surgailis (2009), 
Didier and Pipiras (2011), and others. Regarding o.s.s. processes in general, see Hudson and 
Mason (1982), Laha and Rohatgi (1981), Sato (1991), Maejima and Mason (1994), Maejima 
(1996, 1998), Meerschaert and Scheffler (1999), Section 11 in Meerschaert and Scheffler (2001), 
Chapter 9 in Embrechts and Majima (2002), Becker-Kern and Pap (2008). For related work on 
operator stable measures, see, for instance, Sharpe (1969), Jurek and Mason (1993), Meerschaert 
and Veeh (1993, 1995), Hudson and Mason (1981), among others. 

In particular, Didier and Pipiras (2011) showed that, under the mild assumption 

0<3fi(%)<l, A: = l,...,n, (1.2) 

on the eigenvalues /i^ of the matrix exponent any OFBM Bh admits the so-called integral 
representation in the spectral domain, 

{Bn{t)}tm = { I '-—^{xl^''-'^'^ A^x-}''-^^'^^)B{dx)\ . (1.3) 
Here, x± = max{±x, 0}, 

A = Ai+ iA2 (1.4) 

is a complex- valued matrix with real- valued Ai,A2, A indicates the complex conjugate of A, 
B{x) = Bi{x) + iB2{x) is a complex- valued multivariate Brownian motion satisfying Bi{—x) = 
Bi{x), B2{—x) = —B2{x), and Bi and B2 are independent with induced random measure 
B{dx) satisfying EB(d,x)B{dx)* = dx. Thus, according to (1.3), OFBMs are characterized 
(parametrized) by the matrices H and A. 

In this work, we continue the systematic study of OFBMs started in Didier and Pipiras 
(2011). We now tackle the issues of the symmetry structure of OFBMs and of the non-uniqueness 
(multiplicity) of the exponents which, in our view, are essentially unexplored. Such issues are 
strongly connected. Since the fundamental work of Hudson and Mason (1982), it is well known 
that one given o.s.s. process X may have multiple exponents. More specifically, if we denote the 
set of exponents of X by S{X), we have that 

£{X) = H + T{Gx), (1.5) 
where H is any particular exponent of the process X. Here, 

Gx = {ce GL{n) : {X(i)W = {CX{t)}tm} (1-6) 

is the so-called symmetry group of the process X (where GL{n) is the multiplicative group of 
invertible matrices), and 

TiGx) = \C :C= lim for some {C„} C Gx, < dn ^ o] (1.7) 

L n-i-oo dfi J 

is the tangent space of the symmetry group Gx- By a result for compact groups (e.g., Hoffman 
and Morris (1998), p. 49, or Hudson and Mason (1982), p. 285), it is known that 

Gx = WOqW-^ (1.8) 
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for some positive definite matrix W and some subgroup Oq of the orthogonal group. As a conse- 
quence, the knowledge about (1.5) is subordinated to that about the symmetry group Gx of X. 
For example, the exponent is unique for the process X if and only if the symmetry group Gx is 
finite. 

The description and study of symmetry groups beyond the decomposition (1.8) is a reputedly 
difficult and interesting problem (see, for instance, Billingsley (1966), p. 176, and Jurek and Mason 
(1993), p. 60, both in the context of random vectors; sec also Mccrschacrt and Veeh (1993, 1995)). 
In this paper, we take up and provide some answers for this challenging problem in the context 
of OFBMs. The main goal of this paper is two-fold: to study the symmetry groups of OFBMs 
in as much detail as possible, and based on this, to closely examine (1.5) for OFBMs X = Bh- 
We emphasize again that, to the best of our knowledge, this is the first work where symmetry 
groups are examined for any large class of o.s.s. processes (e.g., for the notion of symmetry groups 
of Markov processes, see Liao (1992) and references therein). Indeed, since its publication, the 
scope of the work of Hudson and Mason (1982) appears to have remained only of general nature, 
the same being true for the main result (1.5). 

The integral representation (1.3) provides a natural and probably the only means to consider 
(almost) the whole class of OFBMs. Section 3 is dedicated to the reinterpretation and explicit 
representation of symmetry-related constructs in terms of the spectral parametrization H, A. 
One of our main results provides a decomposition of the symmetry groups of OFBMs into the 
intersection of (subsets of) centralizers, i.e., sets of matrices that commute with a given matrix. 
For example, in the case of time reversible OFBMs (corresponding to the case when AA* = AA*), 
we show that the symmetry group Gbh is conjugate to 



Here, G(H) denotes the centralizer of a matrix H in the group 0(n) of orthogonal matrices, i.e.. 



and the matrix-valued function Ha, has the frequency x as the argument and is parametrized by 
H and A. Moreover, which is key for many technical results in this paper, we actually express the 
positive definite conjugacy matrix W in (1.8) in terms of the spectral parametrization. This is a 
substantial improvement over previous works on operator self-similarity, where only the existence 
of such conjugacy is obtained, e.g., as in (1.8). 

In view of (1.9) and (1.10), it is clear that the symmetry structure of OFBM is rooted in 
centralizers. The characterization of the commutativity of matrices is a well-studied algebraic 
problem (e.g., MacDuffee (1946), Taussky (1953), Gantmacher (1959), Suprunenko and Tyshke- 
vich (1968), Lax (2007)). We apply the available techniques in a variety of ways to provide a 
detailed study of the symmetry groups and the associated tangent spaces (Sections 4 and 5), as 
well as of the consequences of the non-uniqueness (non-identifiability) of the parametrization for 
integral representations (Section 6). 

Our study of the symmetry structures of OFBMs is carried out from two perspectives: first, 
by looking at the extremal cases, i.e., maximal and minimal symmetry for arbitrary dimension, 
and second, by conveying a full description of all symmetry groups in the lower dimensions n = 2 
and ra = 3. 

Section 4 is dedicated to the first perspective. We completely characterize OFBMs with 
maximal symmetry, i.e., those whose symmetry groups are conjugate to 0{n). We establish 
the general form of their covariance function and of their spectral parametrization. However, 




(1.9) 



G(H) = {O G 0(n) : OH = HO}, 



(1.10) 
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as intuitively clear, maximal symmetry corresponds to a strict subset of the parameter space 
of OFBMs. In view of this, one can naturally ask what the most typical symmetry structure 
for OFBMs is, in a suitable sense. A related question is whether the multiplicity of exponents 
(and, thus, the non-identifiability of the parametrization) is a general phenomenon. Section 4 
contains our answer to both questions, which is, indeed, one of our main results. We prove that, 
in the topological sense, OFBMs with minimal symmetry groups (i.e., {/,—/}) form the largest 
class within all OFBMs. As a consequence, in the same sense, OFBMs generally have unique 
(identifiable) exponents. To establish this result, in our analysis of the centralizers G{Ilx), we 
bypass the need to deal with the major complexity of the eigenspace structure of the function Il-r 
by looking at its behavior at the origin of the Lie group (i.e., as x — > 1), where a great deal of 
information about is available through the celebrated Bakcr-Campbell-Hausdorff' formula. 

Section 5 contains a full description of the symmetry structure of low-dimensional OFBMs, 
namely, for dimensions n = 2 and n = 3. We provide a classification of OFBMs based on their 
symmetry groups. For example, when n = 2, the symmetry group of a general OFBM can be, up 
to a conjugacy, of only one of the following types: 

(i) minimal: {/,—/}; 

(a) trivial: {/, — /, R, —R}, where i? is a reflection matrix; 

(iii) rotational: 5*0(2) (the group of rotation matrices) ; 

(iv) maximal: 0{2). 

Such classification of types for n = 2 stands in contrast with the situation with random vectors, 
for which SO{n) cannot be a symmetry group (Billingslcy (1966)). Nevertheless, we show that 
the latter statement is almost true for OFBMs, since SO{n) cannot be a symmetry group if n > 3. 
In both n = 2 and n = 3, we provide examples of OFBMs in all identified classes, and also discuss 
the structure of the resulting exponent sets £{Bh)- 

In Section 6, wc examine the consequences of non-identifiability for integral representations of 
OFBMs. We show that the multiplicity of the exponents H does not affect the parameter A in 
(1.3) in the sense that the latter can be chosen the same for any of the exponents. Intriguingly, 
this turns out not to be the case for the parameters in the time domain representation of OFBMs, 
and points to one advantage of spectral domain representations. 

In summary, the structure of the paper is as follows. Some preliminary remarks and notation 
can be found in Section 2. Section 3 concerns structural results on the symmetry groups of 
OFBMs. OFBMs with maximal and minimal symmetry groups are studied in Section 4. The 
classification of OFBMs according to their symmetry groups in the lower dimensions n = 2 and 
n = 3 can be found in Section 5. Section 6 contains results on the consequences of the non- 
uniqueness of the parametrization for integral representations. The appendix contains several 
auxiliary facts for the reader's convenience. 

2 Preliminaries 
2.1 Notation 

We shall use throughout the paper the following notation for finite-dimensional operators (matri- 
ces). All with respect to the field M, M{n) or M(7t,,M) is the vector space of all n x n operators 
(endomorphisms), GL{n) or GL{n^M) is the general linear group (invertible operators, or auto- 
morphisms), 0{n) is the orthogonal group of operators O such that OO* = I = 0*0 (i.e., the 
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adjoint operator is the inverse), SO{n) C 0(n) is the special orthogonal group of operators (ro- 
tations) with determinant equal to 1, and so{n) is the vector space of skew-symmetric operators 
(i.e., A* = —A). Similarly, M(m,n, R) is the space of m x n real matrices. The notation will 
indicate the change to the field C. For instance, M(n, C) is the vector space of complex endo- 
morphisms. Whenever it is said that A e M{n) has a complex eigenvalue or eigenspace, one is 
considering the operator embedding M(n) M(n, C). U{n) is the group of unitary matrices, 
i.e., UU* = I = U*U. 5(n,M) is the space of symmetric matrices. We will say that two cndo- 
morphisms A,Be M{n) are conjugate (or similar) when there exists P G GL{n, C) such that 
A = PBP-^. In this case, P is called a conjugacy. The expression diag(Ai, . . . , A„) denotes the 
operator whose matrix expression has the values Ai , . . . , A„ on the diagonal and zeros elsewhere. 
We make no conceptual distinction between characteristic roots and eigenvalues. We also write 
S"'~^ := {v G M" : \\v\\ = 1}; in particular, we denote the complex sphere by S"^"-"^. represents 
a matrix of zeroes of suitable dimension. Whenever necessary, we will specify the dimension of 
the identity matrix by writing /„. Unless otherwise stated, we consider the so-called spectral 
matrix norm ||-||, i.e., \\A\\ is the square root of the largest eigenvalue of A*A. For 
A G M(n, C), we write An — >■ A when \\An — A\\ — )• 0. is the subspace perpendicular to a 
given vector subspace V. For a set of (column) vectors vi, . . . , Vn, A := {vi, . . . , Vn) is the matrix 
whose columns are such vectors. We denote the i-th Euclidean vector by e^, i = 1, . . . , n. 
Throughout the paper, we set 

D = H-^I, (2.1) 

for an operator exponent H. We shall also work with the real part ^{AA*) = ^1^4* -|- A2A2 and 
the imaginary part $J(^^*) = A2AI — A1A2 of AA*. For the real part, in particular, we will use 
the decomposition 

^{AA*) = SbAIS*r = W^, (2.2) 
with an orthogonal Sr, a diagonal Ar and a positive (semi-)definite 

W = SrArS*^. (2.3) 

We shall use the assumption that 

^{AA*) has fuU rank, (2.4) 

in which case Ar in (2.2) has the inverse A^^ . As shown in Didier and Pipiras (2011), the condition 
(2.4) is sufficient (though not necessary) for the integral in (1.3) to be proper and hence to define 
an OFBM. 

All through the paper, we assume n > 2. 

2.2 Remarks on the multiplicity of matrix exponents 

In this section, we make a few remarks to a reader less familiar with the subject of this work. 
It may appear a bit surprising that an o.s.s. process may have multiple exponents, as formalized 
in (1.5). This can be understood from at least two inter-related perspectives: the properties of 
operator (matrix) exponents and the distributional properties of o.s.s. processes. From the first 
perspective, consider for example matrices of the form 

^^={-s o)^"^(2), (2.5) 
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where s G M. Being normal, these matrices can be diagonahzed as Lg = [/2AsC/|, where U2 ^ U (2) 
and As = diag{is, —is). In particular, exp(L2,rfe) = I, k e 7,, since e'^'^'^ = 1. Since Lg and Lg' 
commute for any s, s' G M, this yields 

exp (Lg) = exp (L27rfc) exp (Lg) = exp (L27rfc + Lg), (2.6) 

and shows the potential non-uniqueness of operator exponents stemming from purely operator 
(matrix) properties. Note also that the situation here is quite different from the 1-dimensional 
case: in the latter, the same is possible but only with complex exponents, whereas here the 
matrices L27rfc have purely real entries. 

From the perspective of distributional properties, we can illustrate several ideas through the 
following simple example. The OFBMs in this example will appear again in Section 4 below. 

Example 2.1 (Single parameter OFBM) Consider an OFBM Bh with covariance function 
EBH{t)BH{s)* =: r(t, s) = r/i(t, s)/, where T}i{t,s) is the covariance function of a standard 
univariate FBM with parameter h G (0,1). This process is o.s.s. with exponent H = hi, 
and will be called a single parameter OFBM. Since Bh is Gaussian, O G if find only if 

OT{t, s)0* = r(t, s). In the case of a single parameter OFBM, this is equivalent to OO* = I or, 
since O has an inverse (Bh is assumed proper), 00* = 0*0 = I. In other words, Gbh = 0{n) 
and 

£{Bh) = H + T{0{n)) = H + so{n). 

Thus, a single parameter OFBM has multiple exponents. From another angle, for a given c > 
and L G so{n), we have Llog(c) G so{n) and hence exp(Llog(c)) = G 0{n) = Gbh- Then, 

{BH{ct)}tm = {c'^BHrntm = {c^c^5H(t)}teR = {c^+^5H(t)}teR, 

which also shows that the exponents are not unique. 

For later use, we also note that an equivalent way to define a single parameter OFBM is to 
say that it has the spectral representation 

{BnmteR = {C [ '^^\x\-('^-'^'>B{dx)} (2.7) 
where C is an appropriate normalizing constant and B(dx) is as in (1.3). 
2.3 Basics of matrix commutativity 

We now recap some key facts and results about matrix commutativity that are repeatedly used 
in the paper. To put it shortly, two matrices A,B ^ M{n, C) commute if and only if they share 
a common basis of generalized eigenvectors (see Lax (2007), p. 74). This means that there exists 
a matrix P G GL{n,C) such that we can write A = PJaP~^ and B = PJbP~^, where and 
Jb are in Jordan canonical form. In particular, if A, B are diagonalizable, then they must share 
a basis of eigenvectors. When, for example, A = I, we can interpret that A commutes with any 
B = PJbP~^ e M{n,C) because for (any) P G GL{n,C), A = PIP-^. 

A related issue is that of the characterization of the set of all matrices that commute with a 
given matrix A, the so-called centralizer C{A). In particular, one is often interested in constructing 
the latter based on the Jordan decomposition of A. 

Before enunciating the main theorem on C{A), we look at an example adapted from Gant- 
macher (1959). 
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Example 2.2 Assume the matrix A G M(10, C), with Jordan representation A = PJaP ^, has 
the elementary divisors (i.e., the characteristic polynomials of the Jordan blocks) 



(A - Ai)3, (A - Ai)2, (A - A2)^ (A - A3), (A - A3), (A - A3), 



(2.8) 



where the eigenvalues Ai, A2, A3 are pairwise distinct. Then, C{A) consists of matrices of the form 
X = PXP-^, where 



b a 
c b a 



M(10,C) 3X = 
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(2.9) 



t J 



The blocks on the diagonal correspond to the Jordan blocks of Ja and the empty entries are 
zeroes. 

We now turn to the structure of the blocks for the general case, and (2.9) serves as an illus- 
tration of the latter. We say that a matrix Xa,/3 G M(pa,Pi3, C) has regular lower triangular form 
(e.g., as each of the blocks in (2.9) with letters a through j) if it can be written as 



Xa/3 = 



(^Pc>0), ifpa<P(3, 



(2.10) 



where Tp^ G M{pa,C) is a Toeplitz lower triangular matrix. Also denote by Np^ G M{pa,C) the 
nilpotent matrix 

/O \ 
1 



N, 



Pa 



V 



1 oy 

The next theorem characterizes C{A). The proof can be found in Gantmacher (1959), p. 219 
(see also pp. 220-224). 

Theorem 2.1 Let A G M{n,C), where A = PJaP~^ and Ja is in Jordan canonical form, i.e., 

Ja = diag(AiIpi +Np^,..., A„/p„ + NpJ 

with (not necessarily distinct) eigenvalues Ai, . . . , A„. Then, the general solution to the equation 
AX = XA is given by the formula X = PXjj^P~^ , where Xj^^ is the general solution to the 
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equation JaXj^ = Xj^^Ja- Here, Xj^ can he decomposed into blocks X^p € M{pa,pp,C), a,(3 = 
1, . . . ,u, where 

Y = j ^' if ^13, 

'^^ \ as in (2.10), ifXa = \i3- 

Corollary 2.1 Let A G S{n,R), where A = OTO*, O G 0{n), T = diag(r/i, . . . , r/„). Assume 
r]j^ = . . . = r]ji_ (possibly k = 1) for some subset of eigenvalues of A, and rjj-^ ^ rji for any other 
eigenvalue rji of A. If another matrix M G M(n, M) commutes with A, then M can be represented 
as 

M = {oj^,..., Oj^ , Oil , . . . , Oi^_k )diag(if 11 , K22) (oji , . . . ,Oj^,Oi^, . . . , Oi^_^ )* , 

where diag(ifii, 1^22) is block- diagonal with Ku G M{k,M.), K22 G M(n — A;,R), where 

, . . . , Oji^ are the column vectors in O associated with the eigenvalues rjj-^, . . . , r]j^ , respectively, 
and Oi^,... ,Oi^_i_ are the column vectors in O associated with the eigenvalues r]i^,. ■ ■ ,r]i^_^, re- 
spectively. 

Consequently, spaia.^{ojj^, . . . ,0^,.} is an invariant subspace with respect to M. 

In view of Theorem 2.1 (and Corollary 2.1), it is intuitively clear that, if a matrix F commutes 
with two matrices A and B which exhibit completely different invariant subspaces, then F can 
only be a multiple of the identity. This is accurately stated for the case of symmetric matrices in 
the next lemma, which is used in Section 4. 

Lemma 2.1 Let 

yl, i? G 5^ := {S G 5(n, M) : S has pairwise distinct eigenvalues}. (2-11) 

Assume A and B have no k-dimensional invariant subspaces in common for k = 1, . . . ,n — 1. If 
M G M(n,]R) commutes with both A and B, then M is a scalar matrix. 

The proof of Corollary 2.1 and Lemma 2.1 can be found in Appendix A, together with some 
additional results on matrix commutativity and matrix representations. 



3 Symmetry groups of OFBMs 

Consider an OFBM Bh with the spectral representation (1.3). In this section, we provide some 
structural results on the nature of the symmetry group (see (1.6)). In particular, we explicitly 
express it as an intersection of subsets of centralizers. 

For notational simplicity, denote Gb^ by Gh- Since OFBMs are Gaussian and two Gaussian 
processes with stationary increments have the same law when (and only when) their spectral 
densities are equal a.e., we obtain that 

Gh = {C eGL{n):EBHit)BHis)* = E{CBH{t)){CBH{s))*, s,teR} 
= {C G GL{n) : {x'^A + xZ^A){x-^A + xZ^A)* 

= C{x-^A + xZ^A){x-^A + xZ^A)*C*, xeR} 
= {Ce GL{n) : x-^AA*x-°' = Cx'^ AA*x-^' C* , x > Q} 

= GH,lf^GH,2, (3.1) 
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where 

Gh,i = {C eGL{n) :x'^^{AA*)x-^' =Cx-^^{AA*)x-^''C*, x>0}, (3.2) 
Gh,2 = {C e GL{n) : x-^^{AA*)x-^* = Cx-^'^{AA*)x-^'' C* , x > 0}. (3.3) 

Consider first the set Gh,i- Using the decomposition (2.2) and working under the assumption 
(2.4), we have that 

Gh,i = {C e GL{n) : x'^ SRK\S*jiX-^* = Cx'^ SrK\S*rX-^' C* , x > 0} 

= {C G GL{n) : {A^^ S*rx''Cx-'' SrAr){A^^ S*jix''Cx-'' SrAr)* = I, x>0} 

= {C e GL{n) : A^^ S*rX^Cx-^ SrAr g 0{n), x > 0}. (3.4) 

Taking x = 1 and using the fact that Sr is orthogonal, C G Gh,i necessarily has the form 

C = SRARS*jiOSRA-^S*ji = WOW-^ (3.5) 

with O G 0{n) (see also Remark 3.1 below). Substituting (3.5) back into (3.4), we can now 
express Gh,i as 

Gh,i = W{OeO{n): 

0{W-^x-^^{AA*)x-^*W-^) = {W-^x-^^{AA*)x-^*W-'^)0, x > 0}W-^ 
= Wf]GiU,)W-\ (3.6) 

x>0 

where we use the definition (1.10) of G{Ilx), and 

:= W-^x~^^{AA*)x-^''W-^ = x-^x-^* (3.7) 

with 

M = W-'^DW. (3.8) 

Remark 3.1 A simpler way to write (3.5) and (3.6) would be to replace W = SrArS'^ by SrAr. 
Note that, with our choice, W is positive definite. The relation (3.6) then takes the form (1.8). 

The relation (3.6) describes the first set Gh,i in the intersection (3.1). Instead of describing 
the second set Gh,2 separately, it is more convenient to think of the latter as imposing additional 
conditions on the elements of Gh,i- In this regard, observe first that, for any y > 0, 

Gh,i = y''GH,iy-'', (3.9) 

which simply follows by observing that the condition 

x-^^{AA*)x-^' = Cx-^^{AA*)x-^*C\ x > 0, 

defining the set Gh,i, is equivalent to 

x-^^{AA*)x-^'' = {y^Cy-^)x-^^{AA*)x-^'' {y-^* C*y^'), x > 0. 

Using the relation (3.9), C G Gh,i satisfies the relation 

x-^9(AA*)x-^* = Cx-^'^{AA*)x-'^* C\ X > 0, (3.10) 
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defining the set Gh,2, if and only if C G Gh,i satisfies the same relation (3.10) with x = 1. 
Considering the form (3.5) of C G Gh,i, this imposes additional conditions on the orthogonal 
matrices O. Substituting (3.5) into the relation (3.10) with x = 1, we obtain that 

Q{AA*) = WOW-^Q{AA*)W-^0*W, 

i.e., 

OW-^QiAA*)W-^ = W-^^{AA*)W-^0, 

or 

O e G{Ui), (3.11) 

where 

Ui = W-^Q{AA*)W-\ (3.12) 

By the expressions (3.1), (3.6) and the discussion above, we arrive at the following general result 
on the structure of symmetry groups of OFBMs, and in particular, on the form of the conjugacy 
matrix W. 



Theorem 3.1 Consider an OFBM given by the spectral representation (1-3), and suppose that 
the matrix A satisfies the assumption (2.4)- Then, its symmetry group Gh can he expressed as 

= w^(f| G(n^)nG(n,))vr-\ (3.13) 

where W is defined in (2.3), and Hx and II/ are given in (3.7) and (3.12), respectively. 

The intersection over uncountably many x > in (3.13) can be replaced by a countable 
intersection in a standard way. We have O G f~^x>oG{Iix) if and only if 

Ox-^x-^* = x-^x-^*0, X > 0. (3.14) 
Writing x'^ = Efclo -^''(^ log(^))''/^'' ^^e relation (3.14) is equivalent to 
V V OM'^HM*)'"' ^~ ^"g^^))""' ^~ log(3;))^^ = V V fiin^^O ^~ log(x))'i (- log(x))fe^ 

fel=0fc2=0 ^ ^ fel=0fe2=0 ^ ^ 

or, with ki = k, ki + k2 = m, 

CO rn ^ oc rn ^ 

^ O ^ M^M*)-^ _ (- log(x))- = J2 E MHM*r-'0 _ (- log(x))-. 

m=0 A;=0 ^' m=0 fc=0 ^' 

Equivalently, 

On^'") = n("^)0, m > 1, (3.15) 

where 

rn / X 

Theorem 3.1 can now be reformulated as follows. 



10 



Theorem 3.2 Consider an OFBM given by the spectral representation (1.3), and suppose that 
the matrix A satisfies the assumption (2.4)- Then, its symmetry group Gh can he expressed as 

oo 

Gh = w(^[^ G(n("^))nG(n7))M^-S (3.17) 

m=l 

where W is defined in (2.3), and H^"*) and H/ are given in (3.16) and (3.12), respectively. 

Remark 3.2 Note that the matrix ILx in (3.7) is positive definite. On the other hand, the matrix 
n("*) in (3.16) is symmetric because so are the terms 

('f\M^{M*)'^-^ +( ^ \M'^-^{M*f 

\ fv / \1TI fv / 

defining H^™). However, II^'") is not positive definite in general. For example, with ?R.(AA*) = I 
and normal D, we have 

nM = (£) + £»*)"», (3.18) 

which is not positive definite (not even for m = 1). Note also that 11/ is skew-symmetric, hence 
normal and diagonalizable. 



4 On maximal and minimal symmetry groups 

An operator self-similar process X is said to be of maximal type, or elliptically symmetric, if its 
symmetry group Gx is conjugate to 0(n). At the other extreme, a zero mean (Gaussian) o.s.s. 
process is said to be of minimal type if its symmetry group is {/,—/}. We shall examine here 
these symmetry structures in the case of OFBMs. First, we characterize maximal symmetry in 
terms of the spectral parametrization of OFBMs. Second, we analyze minimal symmetry OFBMs 
through a topological lens. 

4.1 OFBMs of maximal type 

The following theorem is the main result of this subsection. Recall the definition of single param- 
eter OFBMs in Example 2.1. 

Theorem 4.1 Consider an OFBM given by the spectral representation (1.3), and suppose that 
the matrix A satisfies the assumption (2.4). If an OFBM is of maximal type, then it is a single 
parameter OFBM up to a conjugacy by a positive definite matrix. Moreover, this happens if and 
only if 

'^{AA*) = 0, -{D - dI)^{AA*) = ^{AA*){D* - dl), (4.1) 

for some real d. 

Remark 4.1 Conversely, an OFBM which is a single parameter OFBM (up to a positive definite 
conjugacy) is of maximal type (see Example 2.1). We also point out that we have a proof of the 
first claim in Theorem 4.1 which does not make use of spectral representations and dispenses with 
the assumption (2.4). In the proof of Theorem 4.1 we use spectral representations in order to 
illustrate how the main results of Section 3 can be used. 
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Proof: Since the OFBM is of maximal type, the representation Gh = WiO{n)Wi holds, 
where Wi is positive definite. On the other hand, by (3.13), we also have Gh = W20W2^, 
where W2 is positive definite and O is a subgroup of the orthogonal group. Thus, by Lemma A.l, 
O = 0{n). Therefore, in view of Proposition A.l, maximal type occurs if and only if 

= Xa:I, X > 0, Ui = XI, X^, A G M. (4.2) 

Note that 

Ui = XI ^ ^{AA*) = X?R:{AA*) ^ A = o ^{AA*) = 0. (4.3) 

Moreover, 

= X^I, x>0 ^ X^?lt{AA*) = x-^3fi(AA*)x-^*, (4.4) 
which implies that, for any xi,X2 > 0, 

X^,.,,mAA*) = (,xi,X2)-^»(ylyl*)(xiX2)-^* = x^^A,,3ft(^^*)x^^* = X^,X^,^{AA*). 

Hence, under assumption (2.4), Xxix2 = ^xi^x2j xi,X2 > 0. Moreover, the function log(Aexp(.)) 
is additive over M, and it is measurable (since it is continuous). As a consequence, by Theorem 
1.1.8 in Bingham et al. (1987), p. 5, there exists a real d such that log(Aexp(.)) = — 2d(-), i.e., 
Xx = x~^'^. In particular, 

x-^?li{AA*)x-^' = x-'^'^^{AA*). (4.5) 
Relations (4.3) and (4.5) imply that the covariance structure of OFBM can be written as 

r pitx 1 p—isx 1 

EBH{t)BH{sy = / — : \x\-^'^^{AA*)dx. 

In view of the relations (2.2) and (2.7), this shows that Bh is a single parameter OFBM up to a 
conjugacy. 

Finally, note from the above that n^; = Xxl, x > 0, is equivalent to x~^?fi{AA*)x~^* = 
x-^'^^iAA*) or x^-'^^^{AA*) = ^{AA*)x-^^''-'^^^ for x > and some real d. The latter is 
equivalent to {D - dI)^{AA*) = -^{AA*){D* - dl) for some real d. □ 

Remark 4.2 Theorem 6 in Hudson and Mason (1982) shows that every maximal symmetry 
o.s.s. process has an exponent of the form hi, h e R. For the case of OFBMs, the proof of 
Theorem 4.1 retrieves this result (see expression (4.5)). Moreover, it is clear that, for a maximal 
symmetry OFBM Bh (or, as a matter of fact, for any maximal symmetry o.s.s. process) and for 
any H G £(Bh), W~^HW is normal, since W~^(H — hI)W G so(n) (see also Section 5 for further 
results on the structure of exponents for dimensions n = 2 and n = 3). 

4.2 OFBMs of minimal type: the topologically general case 

In view of Theorem 3.1, an OFBM is of minimal type if and only if 

f|G(n,)nG(n,) = {/,-/}, 

x>0 

and, in particular, if 

[]G{Ux) = {I,-I}. (4.6) 

a;>0 
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We shall focus here on the relation (4.6) with the following related goals in mind. 

The first goal is to provide (practical) conditions for (4.6) to hold and, hence, for an OFBM to 
be of minimal type. This is a non-trivial problem. The structure of G{Ilx) depends on both the 
eigenvalues and the invariant subspaces of Tlx, which are arbitrary in principle. Moreover, their 
explicit calculation becomes increasingly difficult with dimension. To shed light on (4.6), we take 
up an idea from Lie group theory: a lot of information about M in the expression = x~^x~^* 
(see (3.7)) is available at the vicinity of the identity in the Lie group, i.e., as x — > 1. The general 
approach we take is to study the behavior of the logarithm of lix through the Baker-Campbell- 
Hausdorff formula, i.e., by looking at the associated Lie algebra. The characterization of the 
behavior of the eigenvectors of will then be retrieved by turning back to the Lie group through 
the exponential map. 

Initially, our conditions for the relation (4.6) to hold are in terms of the matrix M, and not 
directly in terms of H and A. Our second goal in this section is to show that these conditions on 
M yield "most" OFBMs in terms of the parametrization M, and then relate them back to H and 
A. The term "most" is in the topological sense, i.e., except on a meager set. This result should not 
be surprising: if r\x>oG{Jlx) has non-trivial structure, then this imposes extra conditions on M 
(or D, W) as in Section 4.1. Though not surprising, formalizing this fact is not straightforward, 
as shown here. This second goal leads to the main result of this section, which, for the sake of 
clarity, we now briefly describe. In analogy with the assumption (1.2), consider the set 

V = [d€ M(n,M) : < 3?(4) < ^, = 1, . . . , n}, (4.7) 

where di, . . . ,dn denote the charateristic roots of D. Theorem 4.2 below states the existence of a 
set M C M(n, M) such that, for all D and positive definite W such that W~^DW € MnV, the 
OFBM with spectral parametrization D and ^{AA*) := has minimal symmetry. Moreover, 
7W n P is an open set (of parameters), and it is dense in V. Therefore, fl P is a meager 
set. Conversely, every M & M DD gives a minimal symmetry OFBM throTigh an appropriate 
spectral parametrization. In order to provide easy access to the mathematical content of this 
section. Remark 4.6 below contains a short heuristic proof of a weaker version of this claim, i.e., 
that OFBMs have identifiable (unique) exponents for every parametrization except on a meager 
set. 

The rest of this section is dedicated to developing these ideas, as well as the framework behind 
them. Hereinafter, unless otherwise stated, we impose no restrictions on the eigenvalues of M, 
i.e., the expression Hx = x~^x~'^* is taken for any M G M(n,M). Consider the decomposition 
of the latter space into the direct sum 

M(n,M) =5(n,M)eso(n). (4.8) 

Correspondingly, for M G M(n,M), denote 

M = S + L, (4.9) 

where S = (M -|- M*)/2, L = (M — M*)/2 are, respectively, the symmetric and skew-symmetric 
parts of M. 

The next proposition shows that for an appropriately chosen M, the centralizer of the family 
Ila; is minimal. In the proof, the symbol [•, •] denotes the commutator. Since the point x = 1 is a 
singularity in the sense that all the information about M from 11;^ = x~^x~^* is lost at it, the 
idea is to analyze the behavior of Jlx for a; in a close vicinity of 1. The proof of Proposition 4.1 
is based on the idea that, for a matrix parameter M = S + L, S^S^<^ (S(n,]R) and L G so{n), 



13 



the existence of non-trivial solutions for the matrix equations OIlx = ^xO, x > 0, implies the 
coincidence of some invariant subspace of S and L, which, as we will see afterwards, is a very 
special situation in the topological sense. 

Proposition 4.1 Let M = S + L, where S ^ (see expression (2.11)) and L G so{n) do not 
share k-dimensional real invariant subspaces, k = 1, . . . ,n — 1. Then, CI^-^q G{Ilx) = {I, —I}- 

Proof: Note that M + M* = 25, and that 

[M, M*] = MM* - M*M = {S + L){S - L) - {S - L){S + L) = 2[L, S]. 

Since the mapping M i— t- exp(M) is a C°° homeomorphism of some neighborhood of in the Lie 
algebra of GL{n,M) onto some neighborhood U of the identity I in GL{n,M), then its inverse 
function Log is well-defined on U. Therefore, by the Baker-Campbell-Hausdorff formula, for small 
enough log(x) we have 

Log(exp(- log(a;)M) exp(- log(a;)M*)) = - log(a;)(M + M*) + ^[- log(x)M, - log(x)M*] 

+0(log3(x)) = - log(x)(M + M*) + log\x)^[M, M*] + Oilog\x)) 

(see Hausner and Schwartz (1968), p. 63 and pp. 68-69). 

We would like to show that there exists S > such that the symmetric matrices 

-j^Log(x-^x-^*) = (M + M*)-i^[M,M*] + 0(log2(x)), xeB{l,S)\{l}, (4.10) 

do not share A;-dimensional real invariant subspaces with the symmetric matrix M + M*, k = 
l,...,n- 1. 

Without loss of generality, assume that there exists {xj}, — t- 1, such that each associated 
expression in (4.10) shares a /c(i)-dimensional invariant subspace with M + M*. Since the number 
of possible real invariant subspaces of M -|- M* G 5^(n,M) is finite (by Corollary C.2, they are 
generated by real eigenvectors of M -|- M*), by passing to a subsequence if necessary, we can 
assume that each 

-r^^ogixr^xr^") 
\og{xi) ' ' ' 

shares the same invariant subspace with M + M* . Write 

M + M* = 0{2T)0*, 

where T is diagonal and assume without loss of generality that the invariant subspace in question 
is spauRjoi, . . . , Ok} (i.e., the first k columns of O G 0{n)). This implies that we can write 

0(2T)0* - i^§M((OTO*)L - L{OTO*)) + 0{\og'^{xi)) = OJiO*, 
where by Corollary C.l Jj is block-diagonal of the form 

Ji=(^ "o' j1 ) ' -^ii ^ "^22 e M(n -k,R). 

Therefore, 

2T - i^^M(T(0*LO) - {0*LO)T) = Ji + 0{log\xi)). 
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If L2 := 0*L0, then 



TL2-L2T = {Ji-2T)(- 



V log(xi) 



+ 0(log(x,)). 



(4.11) 



Note that the term Jj — 2T is still block-diagonal of the same form as Ji. For notational simplicity, 
we still write Jj — 2T = diag( J|^, J22)- Therefore, the right-hand side of (4.11) must converge as 
Xi 1. Denote the limit by 

Ju 
J22 



J = 



where Ju € S{k,R), J22 £ S{n 



). Denote L2 = {lij)i,j=i, 



TL2-LoT = 



( Zi2(Ai - A2) /i3(Ai - A3) /i4(Ai - A4) 

^23(A2 — A3) /24(A2 — A4) 

h^{\z - A4) 



v 



(4.12) 

,n, T = diag(Ai, . . . , A„). Then 



^ln(Al — A„) ^ 
hn{^2 — An) 
^3n(A3 — A„) 



^n—l,n (An— 1 





An) 



, (4.13) 



where the entries below the main diagonal are equal to the corresponding ones above the main 
diagonal. From expressions (4.13) and (4.12), we conclude that the upper right (non-square) block 
is zero, i.e., 

/ ^l,jt+l(Al — Afe+i) /l,fe+2(Ai — Ajfc+2) • • • ^i,n(Ai — An) ^ 
^2,A;+l(A2 — Afe+i) /2,fe+2(A2 — Afc+2) • • • ^2,n(A2 — An) 



\ h,k+l{^k — Afc+l) h,k+2{^k — Afe+2) 

Therefore, since the A's are pairwise different. 



h,n{^k ~ An) / 



L12 ■= 



( ^l,fc+l ^l,fc+2 
^2,fc+l ^2,fe+2 

V ^fc,fe+l ^fe,fc+2 



^k,n J 



0. 



Therefore, 



L = 



Lu L12 
L21 L22 



o* = 









L22 



o* 



(note that Lu and L22 may contain zeroes off the main diagonal but this is inconsequential). 
Therefore, spanjgjoi, . . . , o^} is also an invariant subspace of L (contradiction). 



So, fix xq G B{1, 8)\{1} such that Log(X[ 



M-M* 



X; 







log(a;o) (M + M* ) + i log^ (xq ) [M, M*] + 



O(log^(a;o)) has no fc-dimensional real invariant subspaces in common with M + M* or, equiva- 
lently, — log(xo)(M + M*), k = I, . . . ,n — I. Moreover, by Lemma B.l, we can assume that, for 
X € 5(1,(5)\{1}, Log(x-^x-^*) G S^{n, 

'xr^"' ) shares no A;- 



We claim that there exists xi G B{l,S)\{l,xo} such that Log(X]^ 



dimensional real invariant subspaces with Log(x^^Xo ^*)- This can be proved by an argument 

similar to the one by contradiction used above in this proof. In fact, assume that there exists a se- 
quence {xi} such that Log(x^^a;^^ ) shares some fc(i)-dimensional real invariant subspaces with 
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Log(xQ ). Since the latter is symmetric, then by CoroUary C.2 its real invariant subspaces 

are generated by k eigenvectors of Log{xQ^ ). Denote a basis of orthonormal eigenvectors of 
the latter by vi, . . . ,Vn- Since Log{xQ^^ x^^ ) G 5^(n,M) and the number of possible subspaces 
of the form span]g{?;j^, . . . ,Vji^} is finite, then we can assume that the shared invariant subspace 
is the same for all i. For notational simplicity, write it as spanR{t;i, . . . ,Vk}; complete this basis 
with orthonormal vectors Vk+i, ■ ■ ■ ,Vn and let P = {vi, . . . ,Vk, Vk+i, . . . ,Vn) € 0(n). Then, by 
Corollary C.l 

(M + M*) + 0(log(x,)) = --^Log(x-^x-^*) = Pdiag(Jii, J^a)^*- 

Thus, limj^oo Jli, ^^^i^oo J22 exist, which implies that span]g{ui, . . . , Vk} is a real invariant sub- 
space of M + M* (contradiction) . 

Consequently, we also have that x^'^^x^'^ and Xq'^Xq'^ share no /c-dimensional real 
invariant subspaces (since the eigenvectors are the same as those of Log{xi^Xi^ ) and 
Log(xQ*^XQ^ ), respectively, and by Corollary C.2, the fc-dimensional real invariant subspaces 
of symmetric matrices are each generated by a set of k real eigenvectors). Thus, by Lemma 2.1, 
the conclusion follows, i.e., the only orthogonal matrices that commute with both Xq'^Xq^ and 
xr^xr^* are ±1. □ 

Bearing in mind Proposition 4.1, we would like to construct a set of matrices M based on 
which we can apply the proposition, and whose topology we can characterize. We now take a 
closer look at an appropriate set of skew-symmetric matrices. 

Definition 4.1 Let oi, . . . , o„ be a real orthonormal basis of M". Let 

Anvar(oi, . . . , o„) = {L G so{n) : 
there exists a subset Oj^, . . . , oji^, k < n, that generates a real invariant subspace of L}. (4.14) 
Example 4.1 If a, b,c,d E M, then 



L :-- 



/ a b \ 

—a c 
-b -c 

d 

V -d 0/ 



G 'Cinvar(ei,e2,e3,e4,e5) 



since the real subspace spanjjjei, 62, 63} is invariant with respect to L (the same being true for 
span]R{e4,e5}). 

Remark 4.3 The representations of the sets >Cinvar may use different arguments (vectors). For 
{vi,V2) = {0} for every orthonormal pair vi,V2 € M . However, this is inconse- 
quential for the developments in this section. 

Proposition 4.2 below establishes the topological properties of the class of "well-behaved" 
exponents M, i.e., those that will eventually be associated with minimal symmetry OFBMs. Its 
proof is based on the next two lemmas. 

Lemma 4.1 (i) The set is an open, dense set in (the relative topology of) 5(n,M). 

(a) For any orthonormal vectors oi, . . . , o„ in R", the set (>Cinvar(oi, . . . , On)y is an open, dense 
set in (the relative topology of) so{n). 
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Proof: (i) Openness stems from the fact that, for {3^} C 5(n, M) such that — >■ Sq, by 
Lemma B.l, the eigenvalues of Sf- converge to those of Sq. Thus, for large enough k, the latter 
are pairwise distinct. 

For denscness, take G {•S^Y- One can obtain a sequence {Sk} C such that Sk — >■ S'o 
simply by appropriately perturbing the eigenvalues of Sq, for a fixed conjugacy O G 0{n) of 
eigenvectors of Sq. 

(ii) By contradiction, fix a real orthonormal basis Oi,...,o„ and assume that (/2invar)'^ • = 
(Anvar(oi, • • • , o„))'^ is not Open. Then there exists L G (Anvar)'^ such that, for some {Ljjjgi^ C 
'^invarj -^i L. Sincc there are finitely many fe-tuples (oj^, . . . , Oj^,), k = 1, . . . ,n — 1, then we 
can extract a (convergent) subsequence {Lii} for which all Lj/'s share one invariant subspace 
spanjjjoj^, . . . ,Oji^} (i.e., k is not a function of i). This means that we can form an orthogonal 
matrix O := {oj^,. . . , Oj^,oj^^^, ojj, where span^joj^^^ , . . . , ojj = (span^loji, . . . , OjJ)^, 
and by Corollary C.l write 

where L\-^ G so{k),L22 G so{n — k). Since \im.i^^ L\i,]imi^^ must exist, then L G Anvar 
(contradiction) . 

As for denseness, once again fix a real orthonormal basis oi, . . . , o„. Take any L G £invar := 
Anvar(oi, . . . , o„) and Write it in a block-diagonal form as L = Odiag(Lii, . . . , Ljj)(9*, where 
Lii, . . . , Ljj are skew-symmetric matrices. Now form the sequence of matrices {Li} by replacing 
all the zero entries above the main diagonal of L with and correspondingly, the zero entries 
below the main diagonal with (this may include entries in the blocks Ln, • • • , Ljj themselves). 
Then, by Corollary C.l, we must have G jO-f^y^j., and ^ L. □ 

We now define a correspondence (set-valued function) that maps the set 5^ into sets of skew- 
symmetric matrices in the classes (4.14). 

Definition 4.2 Let V be the class of all subsets of a set. Define the correspondence (set-valued 
function) 

l:S^-^ V{so{n)), 

S ^ 1{S) = (Anvar(0l, • • • , 0„))^ (4.15) 

where oi , . . . , represent orthonormal eigenvectors of S. 

Remark 4.4 The correspondence /(•) is well-defined. In fact, for S G S^{n,R) and an associated 
basis of eigenvectors oi, . . . , o„, all the possible representations of £invar(oi5 • • • , o„) are of the form 
'Cmvar(±oi , ■ ■ ■ , ±o„). Now uotc that, for any subset Ojj , . . . , Oj^, , 1 < < n, spanj^joj^ , . . . , Oj^ } = 
span]g{±Oj^, . . . , ±Oj^}, i.e., invariant subspaces can equivalently be expressed in either basis. 

In the next lemma, we show that the graph of the correspondence (4.15) is open and dense in 
the set {S^, so{n)). The topology under consideration is that generated by the open rectangles 

Bs{n,R)i^O,£l) X S5o(„)(Lo,£2), £l,£2 > 0, (4.16) 

where 

^5(n,K)(^o,ei) :={Mg M(n,M) : ||M - 5o|| < ei} n 5(n, M), 
Bso(n){Lo,e2) := {M e M{n,R) : ||M-Lo|| < £2} n so(ra) 
and II -11 is the spectral norm. 

Lemma 4.2 Let be the correspondence in Definition 4. 2. Then Graph(Z) := {(S, L) : S G 
S^,L G 1{S)] is open and dense in (5(n, M), so(n)). 
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Proof: Openness is a consequence of the fact that, if Sq G S^, then, as Si — )■ Sq, the eigenvalues 
of Si converge to those of Sq (in the sense of Lemma B.l). Indeed, assume by contradiction that 
there exists {So,Lo) G Graph(/) such that, for some sequence {Si, Li) ^ Graph(Z), 

{Si, Li) — >■ {So,Lo). 

Note that there cannot be a subsequence {Si/} C such that Si' Sq (since this contradicts 
the openness of estabhshed in Lemma 4.1). Thus, we can assume that {SijieN ^ Conse- 
quently, we must have Lj € /3invar(oi5 ■ ■ ■ , o^)) where o|, . . . , o*^ is a basis of real eigenvectors of Si. 
Since Sq G S^, then by Lemma B.2, we can assume that o^^ ^ oi, . . . , ^ o„, where oi, . . . , 
is a basis of real eigenvectors of Sq. Therefore, we can write 

Li = {o\,...,ol)Ki{o\,...,oi,)*, ieN, 

where, by Definition 4.1 and Corollary C.l, possibly by a permutation of the vectors o\,...,d!^ 
the matrix Ki G so{n) can be made block-diagonal (see also Example 4.1). Define 

k* = min{A; = 1, . . . , n — 1 : infinitely many Lj's have a fc-dimensional real invariant subspace}, 

i.e., the minimal non-trivial dimension for invariant subspaces over infinitely many terms of the 

sequence {Li}. Now for the associated sequence of vectors o\,...,ol^, for each i choose one 
subset of indices ji{i), . . . , jk*ii) ^ { 1 ; • • • ? such that o*^ ^.y..., o*^ ^ generates a real invariant 
subspace of Li (there may be more than one choice, but this is inconsequential). Since there 
are at most (j^») such choices, by passing to a subsequence if necessary, we can fix a subset of 
indices ji , . . . , j}.* such that o*-^ , . . . , Oj^^ generates a real invariant subspace of Li for every i in 
this (sub)sequence. Since we can change at will the order of columns in the conjugacy matrix, 
without loss of generality we can assume that ji = 1, . . . = . Therefore, 

Li = 0idiag{L\^,Li2)0*, iGN, 

where L\i G so(fe*), G so{n — k*). Since Oi = {o\, . . . , o*J — >■ (oi, . . . , o„) and Li Lo, then 
the limits limj^oo^in hmj^oo -^22 exist and thus Lq G Anvar(oi, • • • , o„). Therefore, {So,Lq) ^ 
Graph(/) (contradiction). 

We now show denseness. Assume {So,Lo) ^ Graph(Z). We will break up the argument into 
cases. Let oi, . . . , o„ be a basis of real orthonormal eigenvectors of Sq. 

• So G S^: Then Lq G >Cinvar(oi, ■ ■ ■ , On). We can apply the same argument as in the proof of 
Lemma 4.1 and obtain the sequence {So,Li) G Graph(Z), {So, Li) {So,Lq). 

• So ^ S^, Lo G (>Cinvar(oi, . . . , o„))'^: Then generate a sequence Si = ODiO* by appropriately 

perturbing the repeated eigenvalues of So so that Si G and Si — > Sq- Thus, {Si,Lo) G 
Graph(Z) (since Si and So share the eigenvector basis oi, . . . , o„) and {Si, Lo) {So, Lq). 

• So ^ S^, Lo G >Cinvar(oi, . . . , On)'. As in the previous case, generate a sequence Si = ODiO* 
by appropriately perturbing the repeated eigenvalues of Sq so that Si G and Si — t- Sq. 
Without loss of generality, assume that the subset of vectors oi , . . . , Ofc generates a real 
invariant subspace with respect to Lq- Now apply the same argument as in the proof 
of Lemma 4.1 and obtain the sequence {Si, Li) G Graph(Z) (since Si and Sq share the 
eigenvector basis oi, . . . , o^), with {Si, Li) ^ {Sq, Lq). 

□ 
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Remark 4.5 Regarding the last part of the proof of Lemma 4.2 (on denseness), when Sq ^ 
S^{n,'M.) there are infinitely many choices of bases of orthonormal eigenvectors of Sq. Thus, Lq 
is in >Cinvar(oi, . . . , On) Or not depending on the particular basis oi, . . . , o„ chosen. 

The next proposition puts Graph(i) back into the original space M(n, R) in the form of a 
direct sum, rephrases the topological statement of Lemma 4.2, and connects the latter to the 
problem of proving (4.6). 

Proposition 4.2 Let 

M = {M e M{n, M) : M G 5^ e 1{S^)]. (4.17) 

Then, 

(i) M. is an open, dense subset o/M(n,M). Consequently, is a meager set and M is a 
-dimensional C°° manifold in = M(n,M). 

(a) relation (4-6) holds for all M G M.. 

Proof: We first show part (i). Define the norm ||M||^ = HiSH + ||L||, where ||-|| is the spectral 
matrix norm. Expression (4.8) implies that ||-||^ is well-defined. By the equivalence of matrix 
norms, it suffices to show (i) with respect to ||-||0- 

Assume by contradiction that M is not open. Then there is some Mq = Sq + Lq £ M and 
a sequence {M^} C M'^ such that \\Mk — MqH^ 0. However, such convergence holds if and 
only if llS'fc — 5o|| — >■ and ||L/; — Lo|| — >■ 0. Since, for each k, either ^ or ^ l{Sk), then 
this contradicts the openness of Graph(/) (Lemma 4.2). Denseness can be addressed in a similar 
fashion, and the geometric statement is immediate. 

Part (ii) is a consequence of Proposition 4.1. □ 

Example 4.2 To construct an example of M G A^, we turn to the case when n = 3. Take S = 
diag(si, S2, S3), where the (real) eigenvalues are pairwise different. Then, take any L G so(3)\{0} 
not having a real invariant subspace generated by a 1 or 2-tuple of Euclidean canonical vectors 
ei) 62, 63. In other words, we cannot take a matrix L of one of the forms 

OaO\ / b \ /0 0\ 
-a 00 , 00 , c , 

00/ \-6 00/ \0 -c J 

where a, 6, c G M. Now set M = 5 + L. 

In order to make the claim about the general minimality of the symmetry groups of OFBMs, 
we need to restrict the parameter space, as in (1.2). For this purpose, we consider the set V in 
(4.7). The following is the main result of this section. It shows that, except possibly when the 
parametrization is taken on a meager set, OFBMs are of minimal symmetry. 

Theorem 4.2 For all D,W & M(n,M), where W is positive definite, and such that 

W-^DW eMnv, 

the associated OFBM with spectral parametrization D and^{AA*) := W'^ has minimal symmetry. 
The set M.r\V is open, and, in particular, it is an n'^ -dimensional C°° manifold in = M{n, M). 
Moreover, it is also a dense subset ofV (in the relative topology ofV). As a consequence, M'^CiV 
is a meager set. 

Conversely, every M £ M. CiV gives rise to a minimal symmetry OFBM through the spectral 
parametrization D := M, W := I. 
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Proof: By the convergence of eigenvalues ensured by Lemma B.l, !> is an open set. Therefore, 
by Proposition 4.2, fl P is also an open set. The geometric statement is straightforward. 
Furthermore, since A4 is dense in M(n,M), then M. d V must also be dense in the relative 
topology of the open set V. 

The converse is an immediate consequence of Proposition 4.2. □ 

Remark 4.6 A simple heuristic argument may shed light on a weaker version of the claim of 
Theorem 4.2, i.e., with identifiability (uniqueness) of the matrix exponent H in place of the 
minimality of the symmetry group. 

Consider a matrix parameter S + L = M = W~^DW such that S G and L G so{n). For 
X ^ 1, but close to 1, the Baker-Campbell-Hausdorff formula gives 

Logix-^x-^*) = - \og{x){M + M*) + 0(log2(a;)). 

Since 

log(x) 

then the term 0(log^(x)) is only a slight perturbation of the matrix — log(x)(M + M*). Assume 
for simplicity that 0{log^{x)) = 0. Then = eL°g(^"*'^"'^*) = e-^og{x){M+M*) jg symmetric 
with pairwise different eigenvalues. We thus obtain 

G(n,) = p,(±i,±i,...,±i)p,*, 

where Px G 0{n) is a matrix of eigenvectors of n^;. In particular, G{Ilx) is finite and thus the 
matrix exponent of the associated OFBM is identifiable. 

Moreover, the set is open and dense in S{n, M), thus implying that © so{n) is open and 
dense in the subset of M(n,M) whose eigenvalues have real parts between —1/2 and 1/2. 



5 Classification in dimensions n = 2 and n = 3 

Theorems 3.1 and 3.2 describe the general structure of symmetry groups of OFBMs. The cases 
of maximal and minimal symmetry groups were studied in Section 4. In this section, we are 
interested in identifying all the possible "intermediate" symmetry groups. We shall describe their 
structure in dimensions n = 2 and n = 3, and make some comments about higher dimensions. 



5.1 Dimension n = 2 

When n = 2, the contribution of the term G{Ili) in Theorems 3.1 and 3.2 can be easily described, 
as the next two results show. 

Lemma 5.1 IfUi 7^ 0, then G{Ui) = 50(2). 

Proof: Since 11/ is skew-symmetric, we have 

Thus, IIj/A is a rotation matrix, and thus G{Ilj) = S0(2). □ 
Theorems 3.1 and 3.2 can now be reformulated as follows. 
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Corollary 5.1 For n = 2, under the assumptions and notation of Theorems 3.1 and 3.2, we have 

r - / n,>oG(n,)n50(2), if / o 1 i 

" I n.>oG(n,), if5>(AA*) = ^^-'^ 

_ ^ / n^>iG(n(™))n50(2), if9(^A*)/ol ^ 

" n^>iG(nM), if = / ^ • ^^-^^ 

Next, we study the possible structures of the groups G{U) when 11 is symmetric (and hence 
potentially positive definite, as the matrix in (5.1)). Let 7ri,7r2 be the two real eigenvalues of 
n. Two cases need to be considered: 



Case 2.1: vri = 7r2, 
Case 2.2: tti ^ ■n2- 



In Case 2.1, 11 = tti/ and hence 
In Case 2.2, we can write 



(5.3) 

G(n)=0(2). (5.4) 



where the columns of the orthogonal matrix S = {pi P2) consist of the orthonormal eigenvectors 
pi,P2 of n. By Theorem 2.1, B G 0(2) commutes with such 11 if and only if B = SOS* where G 
is a diagonal matrix such that G^ = I (G^ = / is a consequence of the fact that B G 0(2)), or 

We thus have 

G(n) = {i.-i.s(^l \ )^*} 

= {/,-/, Ref(pi),Ref(p2)}, (5.5) 

where Ref(p) indicates a reflection around the axis spanned by a vector p. The expressions (5.4) 
and (5.5) provide the only possible structures for G(n). Together with Corollary 5.1, this leads 
to the following result. 

Theorem 5.1 Consider an OFBM given by the spectral representation (1.3), and suppose that 
the matrix A satisfies the assumption (2.4). Then, its symmetry group Gh is conjugate to one of 
the following: 

(2. a) minimal type: {I, —I}; 

(2.6) trivial type: {I, —I, Ref(pi), Ref(p2)} for a pair of orthogonal pi,P2; 
(2.c) rotational type: 5(9(2); 
(2.d) maximal type: 0{2). 

All the types of subgroups described in Theorem 5.1 are non-empty, as we show next. Since 
OFBMs of maximal and minimal types were studied in general dimension n in Section 4, we now 
provide examples of OFBMs of only the two remaining types for dimension n = 2. 
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Example 5.1 (Rotational type) Consider an OFBM with parameters 

D = dl, V2A1 G 50(2)\{/, -/}, V2A2 G 0(2)\50(2), (5.6) 

where d is real. Then, we have = x-^'^I and G(n^) = 0(2). Since 9(^^*) ^ 0, Corollary 5.1 
yields that Gh = S0{2). 

Example 5.2 (Trivial type) Consider an OFBM with parameters 

where di ^ ^2 are real. Then, 

AA* = ( I'^jl^ 2 ) = ^{AA*), Q{AA*) = 0, 



and 



_ , x-^'^i 



implying that, for x / 1, 

G(n,) = <{/,-/, 

Corollary 5.1 then yields 



1 \ / -1 
-1 / ' I 1 



Only OFBMs of rotational and maximal types have multiple exponents. Moreover, in view of 
(1.5), in both cases we have 

£{Bh) = H + Wso{2)W-^, (5.9) 

where H is any exponent of the OFBM Bh- This relation can be further refined, as the following 
proposition shows. For this purpose, we need to consider a so-called commuting exponent Hq G 
£{Bh), i.e., an exponent Hq such that 

HqC = CHo (5.10) 
for all C G Gh- The existence of this useful exponent is ensured by Lemma 2 of Maejima (1998). 

Proposition 5.1 Consider an OFBM given by the spectral representation (1-3), and suppose 
that the matrix A satisfies the assumption (2.4) ■ If S{Bh) is not unique, then the commuting 
exponents are of the form 

Ho = WU2diag{h,h)U2W-^, (5.11) 

where 

U. = ^(' \) (5.12) 



2 

and h E C In particular, 

£{Bh) = W{U2dia.g{h,h)U2 + so(2))VF-\ (5.13) 
H = ^{h)I G £{Bh) and W'^HW is normal for any H G £{Bh). 
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Proof: If £{Bh) is not unique, then by Theorem 5.1, Hq commutes with WSO{2)W~^. In 
particular, Hq commutes with WOW~^ for O G S'0(2)\{/, — /}. Since such O is diagonahzable 
with two complex conjugate eigenvalues, the eigenvectors of WOW~^ are also eigenvectors of 
Hq. Thus, Hq can be written as VF[/2diag(/ii, /i2)f/2 Therefore, since hi, h2 are also the 
eigenvalues of ?72diag(/ii, /i2)i7|, which must only have real entries, a simple calculation shows 
that hi = h2, and thus (5.11) holds. This also yields (5.13). 

For H G £{Bh), (5.13) implies that W~^HW is normal. In particular, we may choose the 
exponent H = Hq + WL_c^(^,^-^W~'^ = ^{h)I, where Lg is defined in (2.5). □ 

Remark 5.1 In the case of OFBMs of rotational type, which have multiple exponents, every 
exponent is a commuting exponent (compare with Meerschaert and Veeh (1993), p. 721, for the 
case of operator stable measures) . 

Remark 5.2 For general proper Gaussian processes, one can define symmetry sets (groups) in the 
same way as for o.s.s. processes, and, in particular, show that they are also compact subgroups of 
GL{n, R). By applying the argument of the proof of Theorem 4.5.3 in Didier (2007), which is based 
on general commutativity results, instead of spectral filters, one can show that the classification 
provided by Theorem 5.1 actually holds for the wide class of proper bivariate Gaussian processes. 

5.2 Dimension n — 3 

We will make use of the partition of 0(3) into the following subsets: 

SO{3) = {/} U Rote U Rot^, 0(3)\50(3) = {-/} U Refg U Refo, 
where for a vector p, 

Rote ■= IJ R'Ote(p), Rot^ := IJ Rot^{p), 

pes'"-! pg5"-i 

Refe:= IJ Refeip), Refo := |J Refo(p), 

and 

RotQ{p) = {rotations about the axis span]g(p) by an angle not equal to tt}, 
Rot,r(p) = {rotation about the axis span]g(p) by an angle equal to vr}, 

Ref0(p) = {rotations about the axis span]g(p) by an angle not equal to vr, combined with the 

reflection in the plane through the origin which is perpendicular to the axis}, 
Refo(]3) = {reflection in a plane through the origin, where the plane is perpendicular to p}. 

(5.14) 

From a matrix perspective, for some p G S*""^, 

cos 9 — sin \ 
Rote{p)=\ sine cos9 , 6* G (0, 27r)\{7r}, Rot^(p) 
1/ 

COS — sin ^ 

sin 6 COS 6 
0-1 

where = indicates conjugacies by orthogonal matrices. 



1 

1 











-1 











1 










(" 


1 





V 





-1 



23 



Remark 5.3 The subscript 9 in Rotg or Refg only indicates that the angle in question is not 
or TT. Here, 6 does not refer to a specific angle. Indeed, even in the case of a fixed p, Rot5)(p) and 
Ref0(p) are classes of matrices. Also, in the expression Refo we use the subscript to indicate 
that there is no rotation before reflection through the plane in question. 

We first describe the possible structures of G(Ji) for symmetric matrices 11 (such as the 
matrices lix, x > 0, in (3.13)). Let 7ri,7r2,7r3 be the three real eigenvalues of 11. Three cases need 
to be considered, namely. 

Case 3.1: tti = 7r2 = vrs. 

Case 3.2: tti = 7r2 7^ tts, (5.15) 
Case 3.3: tt^ 7^ tTj, i 7^ j, i,j = 1, 2, 3. 

The next proposition gives the form of G{Il) in all the above cases. 

Proposition 5.2 Let 11 G 5(3, M). Denote its eigenvectors by pi, i = 1,2,3, where S = 
{Pi P2 Ps) G 0(3). Then, 

(i) in Case 3.1 in (5.15), 

G{Ii) = 0(3); (5.16) 

(a) in Case 3.2 in (5.15), 

G(n) = {I, -1} U (Rot0(p3) U Reieipz)) U (Rot^(p3) U Refo(p3)) 

U U (Rot,((?)URefo((z)); (5.17) 

?espan,j,{pi,p2} 

(Hi) in Case 3.3 in (5.15), 

G(n) = {/, -/,Refo(m),Refo(p2),Refo(P3),Rot^(m),Rot^(P2),Rot^(p3)}- (5.18) 

Proof: (z) is immediate, so we turn to (m). In this case, we can write 11 = S'diag(7ri, tti, 7r3)S'*. 
By Theorem 2.1, B commutes with such 11 if and only if 

(cii C12 \ 
C21 C22 S*, (5.19) 
c? / 

where C = (c.jj) jj=i^2 and d are arbitrary. If we are only interested in orthogonal matrices, this 
gives C G 0(2), and d = ±1, which corresponds to the subgroup (5.17). Indeed, the matrices 
R'Ot0(p3) and Rot7r(p3) in (5.17) account for rotations C G 50(2) and d = 1 in (5.19), Ref5)(p3) 
and Refo(^33) in (5.17) account for rotations C G 50(2) and d = — 1 in (5.19), and Rot7r(q') and 
Refo((7), q G SY>aii^{pi,p2\, account for reflections C G 0(2)\S'0(2) and d = ±1 in (5.19). 

Regarding {iii), we can write 11 = S'diag(7ri, 7r2, vra)^'*. By Theorem 2.1, B G 0(3) commutes 
with such n if and only \{ B = S'diag(ibl, ±1, ±1)5"*. We thus have 

G(n) = {I, -/, 5diag(-l, 1, 1)5*, 5diag(l, -1, 1)5*, 5diag(l, 1, -1)5*, 
5diag(l, -1, -1)5*, 5diag(-l, 1, -1)5*, 5diag(-l, -1, 1)5*}, 

as stated. □ 
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The expressions (5.16), (5.17) and (5.18) describe the only possible structures for G(n), all of 
them, as shown below, being symmetry groups of some non-empty subclass of OFBMs. However, 
new symmetry groups may arise when one considers intersections of G{Ilx) for different values of 
X, and also with G{Ili). In order to provide a full description of symmetry groups of OFBMs in 
dimension n = 3, we first consider the case of time reversible OFBMs, before turning to the general 
case. As shown in Didier and Pipiras (2011), time reversibility corresponds to the assumption 
that 

'^{AA*) = 0. (5.20) 
Under (5.20), the presence of G(n/) in (3.13) and (3.17) can be ignored. 

Theorem 5.2 Consider an OFBM given by the spectral representation (1-3), and suppose that the 
matrix A satisfies the assumptions (2.4) and (5.20). Then, its symmetry group Gh is conjugate 
by a positive definite matrix W to one of the following: 

(3. a) minimal type: {I, —I}; 

(3.6) for some vector p, 

{I,-/,Refo(p),Rot^(p)}; 

(3.c) for some orthogonal pi,p2,P3, 

{I, -I, Refo(pi), Refo(p2), Refo(p3), Rot^(pi), Rot7r(p2), Rot^(p3)}; 

(3.d) for some orthogonal pi,p2,P3, 

{I, -1} U (Rote(p3) U Refeips)) U (Rot^(p3) U Refo(p3)) U (J (Rot^(g) U Refo(g)); 

?espan,j{pi,p2} 

(3.e) maximal type: 0(3). 

Proof: Recall that, under (5.20), the symmetry group Gh is conjugate to r\x>oG{Ilx). By 
Proposition 5.2, G{Ilx) can only be of the forms (5.16), (5.17) and (5.18). The proof is now split 
into the following cases. 

Case 1: for some x > 0, G{Ilx) has the form (5.18). Since the intersection with some other 
G{Ilx') can only reduce the group, r]x>oG(n.x) can be only of types (3. a), (3.6) or (3.c), where 
(3.6) is a consequence of intersecting (5.18) with some appropriate (5.18). 

Case 2: all G(Jlx) are the same and have the form (5.17). This gives type (3.(1). 

Case 3: there are xi ^ X2 and two different G(nxi) and G(nx2) such that both have the 
form (5.17). Let ri,r2,r3 and Pi,P2,P3 be the corresponding orthonormal vectors in (5.17). For 
G{Ilxi) and G{Ilx2) to be different, the corresponding axis span]g{r3} and spauj^jps} have to be 
different. We break up the proof into subcases. In all of them, based on necessary conditions for 
commutativity, we will construct a set C of "candidate" matrices such that C ^ G{Ilxi) H G{Ilx2)- 
As we will see, in all subcases, C has one of the forms (3. a), (3.6) or (3.c). Since all the possible 
symmetry groups which are subgroups of (3.c) (the most encompassing of the three) are of the 
forms (3. a), (3.6) or (3.c), then nx>oG{Ilx) must be of one of the latter forms, which completes 
the proof. 

Thus, the subcases and the respective sets C are as follows. 
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(3.i) p3 £ spaiiglri, r3}\(spanK{ri}UspanK{r3}) (or with r2 in place of ri): then, by Theorem 2.1 
(c.f. the proof of Proposition 5.2, case (m)), both spanj^jrs} and spanjjjps} are real invariant 
subspaces of any M G G{Ilxi) H G{Ilx2)- Consequently, since and are not orthogonal, 
and two real eigenvectors of orthogonal matrices associated with different eigenvalues 1 
and —1 must be orthogonal, for M we must have that spanR{r3,p3} = span^jri, r3} is a 
two-dimensional real (proper or not) subspace of the eigenspace associated with either 1 
or —1. Therefore, the remaining eigenvalue of M must be real, and by Lemma A. 2 there 
exists an associated eigenvector which is orthogonal to spauj^jri, rs}. In particular, r2 is an 
eigenvector. Therefore, we obtain the set of candidate matrices 

C = {±I,±(ri,r2,r3)diag(l,-l,l)(ri,r2,r3)*}. 

(S.ii) p3 G spanjgjri, r2}: then for M G G{Ilxi) H G{Ilx2), by the same argument as in (3.i), 
both p3 _L r3 are (real) eigenvectors of M associated with real eigenvalues. Thus, the third 
eigenvalue of M is also real and, by Lemma A. 2, {s])an^{p3,r3})^ C spanjjjri, r2} contains 
a unit norm real eigenvector of M, which we can denote by q. Thus, we can set 

C = {M:M = (g,p3,r3)diag(±l,±l,±l)(g,p3,r-3)*}. 

(3.m) ps ^ Uj^jspanjjjrj, rj}: then, by the same argument as in (S.i), spanjgjrs} and span]g{p3} 
are both real invariant subspaces for the solutions. However, since they are not orthogonal, 
then, for each M G G(n.xi) H G(n.x2)i we have that spayi^{p3,r3} is a two-dimensional real 
(proper or not) subspace of the eigenspace associated with either 1 or —1. Thus, by Lemma 
A.2, there exists q G spaujjjps, rs}-*-, ||g|| = 1, such that we can write 

C = {M : M = (9,p3,r3)diag(±l,±/2)(g,P3,r3)"^}. 

Note that the matrices in C are, indeed, orthogonal. In fact, set vi = q and = r^. Now let 
'^ri,r2 be the operator for the projection on spauj^jri, r2}, and set V2 = ir~^^~^T^TTT ■ Then 
V2 7^ 0, because p^ and are not coUinear. Moreover, we have 

^^2'?) = Ju Kr^\\i^ri,r2{P3),<l) = JZ ^7-^(^3,7^^1,^2(9)) = JZ ^7-^(^3, ?) = 0, 

lFri,r2(P3j|| lFri,r2(P3j|| iFn.ra (P3 j || 

where we used the self-adjointness of 7rri,r2 (') the fact that q G spauj^jri, r2}. Therefore, 
we have that V2 G spanjgjps, i;3}. Moreover, since V2 _L rz, then vi,V2,V3 are linearly 
independent and thus form an orthonormal basis. Then we can write 

C = {M:M=(q, t^^^^^^^, n) diag(±l, ±12) (q, /H^' *}• 
^ lFri,r-2(P3)|| ^ IKri,r2(P3)|| 

Case 4: For all x > 0, G{Ilx) has the form (5.16). This gives type (3.e). □ 

We now provide examples of OFBMs of the types (3.6), (3.c) and (3.d), thereby showing that 
all the types described in Theorem 5.2 are non-empty. 

Example 5.3 (Type (3.6)) Consider the OFBM with spectral representation parameters A := I 
and 

/ d \ 
D = i 1 d \ . 
\0 d J 
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By Theorem 3.1, we may assume that the positive definite conjugacy associated with G{Bh) is 
W = I. Observe that 

/ 1 -log(x) 

= x~^x-^* = x-^'^ - log(x) log^(a;) + 1 
V 1 

Due to the block-diagonal shape of x~^x~^* , it suffices to focus on its 2 x 2 upper left block. 
Consider x = e~^,e. The associated 2x2 blocks, i.e., 

1 1 \ / 1 -1 

1 2 y ' \ -1 2 

have pairwise different eigenvalues. Moreover, they do not share eigenvectors, since otherwise 
they would commute. As a consequence, by Proposition 5.2, they are of the form (3.6) with 

p = (o,o,iy. 

Example 5.4 (Type (3.c)) Consider the OFBM with spectral representation parameters 

D = diag(di, 6^2,^3 ), ^ = diag(ai,a2,a3), (5-21) 

where di / dj, i + j. Then, AA* = diag(|aip, |a2p, laap) = ^{AA*), Q{AA*) = and = 
diag(x-2<^i , , ) . This yields 

Gh = {/,-/, Refo(ei),Refo(e2),Refo(e3),Rot^(ei),Rot^(e2),Rot^(e3)} (5.22) 

for the Euclidean vectors Cj, i = 1,2,3. 

Example 5.5 (Type (3.cZ)) Consider the OFBM with spectral representation parameters 

D = disig{di,di,d3), ^ = diag(ai, 02, 03), (5.23) 

where di / ^3. Then, AA* = diag(|ai|Ma2|Ma3p) = ^{AA*), Q{AA*) = and = 
diag(x~^'^i , x"^*^! , a;~^* ) . This yields 

Gh = {I, -1} U (Rote(e3) U Refe(e3)) U (Rot^(e3) U Refo(e3)) U [J (Rot^(g) U Refo(g)). 

gespan{ei,e2} 

(5.24) 

We now extend Theorem 5.2 to the general case of OFBMs which are not necessarily time 
reversible, i.e., we drop the assumption (5.20). From the perspective of the structural result 
provided by Theorem 3.1, the lack of time reversibility manifests itself as an additional constraint 
which may reduce the symmetry group, and even generate a new type, as seen in the next theorem. 

Theorem 5.3 Consider an OFBM given by the spectral representation (1.3), and suppose that 
the matrix A satisfies the assumption (2.4)- Then, its symmetry group Gh is conjugate by a 
positive definite matrix W to the ones described in Theorem 5.2, plus the following: 

(3./) for some vector p, 

{I, -7, Refo(p), Rot^(p), Refe(p), Koie(jp)}. 
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Proof: If 11/ = 0, then G{Ili) = 0{n). So, assume II/ 7^ 0. By the same argument as in 
Theorem 5.2, Case 3, intersecting G(U.j) with any of the subgroups (3. a), (3.6) and (3.c) imphes 
that, eventually, the resulting symmetry group must be of one the forms (3. a), (3.6) or (3.c). 
Therefore, we may only look into the intersection of Gijii) with subgroups of the form (3.(i). 
Assume that the latter are expressed with respect to an orthonormal basis ri,r2,r'3. 

Since 11/ G so(3), then there exists Sj := (pi,P2,P3) G 0(3) such that 11/ = S'/diag(Ls, 0)51, 
where Lg ^ has the form (2.5). Therefore, by Theorem 2.1, we have 

G(Ui) = Sidiag{SO{2),±l)S*i. (5.25) 

Note that, for a matrix of the form Odiag(50(2), ±1)0* = 0^/diag(e*^ 6"*^ ±l)C/*0*, 6 G 
(0, 27r)\{7r}, where U = diag(C/2,l) and U2 is as in (5.12), only the eigenvalue 1 (or — 1) 
is associated with a purely real eigenvector. Thus, we can break up the rest of the proof 
into two cases. Assume that spauj^jps} = spanjgjrs}. Then wc can write G{Ili) = 
(ri,r2,r3)diag(50(2),±l)(ri,r2,r3)*. Thus, G(n/) n (n, r2, r3)diag(0(2), ±l)(ri, rs, rs)* = 
G{Ilj), which gives (3./). Alternatively, assume span]g{p3} / span]g{r3}. In this case, one 
can argue exactly as in the proof of Theorem 5.2, Case 3, {i)-{ni), to obtain the same sets of 
candidates C. Thus, the claim follows. □ 

Example 5.6 (Type (3./)) Analogously to Example 5.1, consider an OFBM with parameters 

D = dl, ^{AA*) = I, '^{AA*) = diag(L, 0), L e so(2)\{0}, 
where d is real. Then, G{Ilx) = 0(3) and G{Ili) is as in (5.25). 

Theorems 5.2 and 5.3 stand in contrast with Theorem 5.1 in that they show the much greater 
wealth of possible symmetry groups in dimension 3 as compared to dimension 2. In a certain sense, 
this enhances the claim of Theorem 4.2 in that, notwithstanding the increasing complexity of the 
possible symmetry structures as dimension increases, minimal type symmetry groups remain the 
topologically general case for any dimension. 

We now provide the tangent spaces and exponent sets for each symmetry group with non- 
trivial tangent space. The proof is along the lines of that for Proposition 5.1. 

Proposition 5.3 Under the assumptions of Theorem 5.3, for the symmetry groups associated 
with non-trivial tangent spaces, the tangent spaces, commuting exponents Hq and sets of exponents 
have the form: 

(3.d) for some orthonormal p\, p2, P3, and the associated matrix S := (pi, p2, Ps), 

T{Gh) = WSdi&giso{2),0)S*W-\ Hq = WSUdiag{hi,hi,h2)U*S*W-\ 

£{Bh) = WSU{diag{hi,hi,h2) + diag{so{2),0))U* S*W-^, 
where U = diag([/2) 1) (md U2 is as in (5.12), and hi G C, /i2 G IK; 
(3.e) T{Gh) = T{SO{3)) = Wso{3)W-^, Hq = HqI, S{Bh) = HqI + Wso{?>)W-^; 
(3.f) the same as for (3.d). 
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Proof: For type (3.d), just note that T{Gh) = T{Rote{p3)) = WSdiag{so{2),0)S*W-^, from 
which Hq and S{Bh) promptly follow. The same argument holds for type (3./). 
The case of type (S.e) is straightforward, since T{Gh) = T{SO{3)). □ 

Remark 5.4 In general dimension n, there are no additional difficulties in describing the struc- 
ture of groups G(n) for a fixed symmetric matrix 11. Equivalently, one can generalize Proposition 
5.2 to the context of dimension n without much effort. Nevertheless, it is cumbersome to describe 
the structure of intersections G(ni) H G(n.2), which is needed for the full characterization of sym- 
metry groups Gh as in (3.13) and (3.17). At this point, a full description of symmetry groups in 
general dimension n is an open question. 

Remark 5.5 The classification given in Theorem 5.1 stands in contrast with the fact that SO{n) 
cannot be a symmetry group for ]R"-valued random vectors (Billingsley (1966)). In particular, 
5*0(2) is not a maximal element of its equivalence class of subgroups in the sense of Meerschaert 
and Veeh (1995), p. 2 (not to be confused with the symmetry group of maximal type in Theorem 
5.1). However, it turns out that Billingsley's result is actually almost true for OFBMs, and 
more generally, proper zero mean Gaussian processes. In other words, for the latter class of 
processes, SO{n) can only be a symmetry group when n = 2 (cf. Theorem 5.3). Indeed, without 
loss of generality, assume W = I. Then, it suffices to show that SO{n) C G{X) implies that 
0(n) = G{X) when n > 3. However, the latter equivalence is a consequence of Proposition A.l 
in the appendix. 

6 On integral representations of OFBMs with multiple exponents 

In this section, we show that when an OFBM has multiple exponents, the matrix A in (1.3) can 
be chosen the same, no matter what matrix exponent is used in the parametrization. We also 
show that, by contrast, such invariance of the parametrization does not hold for the so-called time 
domain representation of OFBM. 

We first consider the latter point. Under (1.2) and 3ft(/i) / 1/2 for any eigenvalue h of H, the 
OFBM {BH{t)}teR also admits an integral representation in the time domain, i.e., 

{BHmteR = I / m - u)"-'^' - (-«);~^')M+ + at - uf-'^' - i-ut-'^')M_)Bidu)} , 

(6.1) 

where M_|_, M_ G M(n, M), and {B(u)}ueR is a vector-valued process consisting of independent 
Brownian motions and such that EB{du)B{du)* = du (Didier and Pipiras (2011)). The following 
example shows that, in general, the matrix parameters M+, M_ cannot be chosen independently 
of the exponent. 

Example 6.1 Consider a bivariate OFBM Bh with the time domain representation (6.1), where 

D = dl, d e (-l/2,l/2)\{0} (or H = hi, h G (0, l)\{l/2}), M+ = O e SO{2) and M_ = /. 
Since rotation matrices commute, it follows directly from (6.1) that 

SO(2) C Gh. (6.2) 

The relation (6.2) implies that r(5'0(2)) = so{2) C T{Gh)- Hence, in view of (1.5), 

H + Lc, c G M, (6.3) 
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axe the exponents of the OFBM Bh, where Lc G so(2) is given in (2.5). Thus, the OFBM Bh 
has the time domain representation 

{BnmteR = { / {{{t - «)?+^^ - (-u)^+^=)M+ + {{t - «)^+^^ - {-u)^+'^^)M_)B{du)} 

(6.4) 

where M+ = M+(c), M_ = M_(c). We want to show that one cannot generally take the original 
parameters M+ = O, M_ = / in the representation (6.4). 

Arguing by contradiction, suppose that = O, M_ = I in (6.4) lead to the same OFBM 
for any c G M. In the spectral domain, these processes have the representation 



/ 

Jr 



where B(dx) is as in (1.3), and 



pitx _ 1 _ 

ixl^-^^A + xZ^-^'A)B{dx), (6.5) 



v27r 



(see Theorem 3.2 and its proof in Didier and Pipiras (2011)). Note that x~^'' commutes with A, 
A and hence {—Lc) can be removed from the exponents of x_ in (6.5). Then, if (6.5) is the 
same process for all c G M, the matrix 

(27:) AA* = T{dl + Lc + I) (e-'<'^^+^^y^O + e'^^'^^+^^^/^I) ■ 

.(Q*gm(d7+L*)/2 ^ ie-^<di+Lt)/2-^-p{dI + Lc + I)* 
= r{dI + L„ + I)r{dI + Lc + iyie-''"^^0 + e''"^^0* + e-'''^-+e'''^^) (6.6) 

does not depend on c. Note that 

T{dl + Lc + I)r{dl + Lc + 1)* = [/2diag(|r(d + ic+ \T{d + ic+ l)f)U2 = \r{d + ic+ 

where U2 is as in (5.12) and T{d + ic + 1) is the univariate Gamma function evaluated at d + 
ic + 1 G C. Writi ng O = ?7 2diag(e^^, e~'^)[/|> for some (3 G (0, 27r)\{7r}, the matrix (6.6) becomes 
C/2diag(/(d, c, P), f{d, c, ^))C/|, where 

fid, c, P) := \T{d + ic + 1)|2 (^e-'^^^e'^ + e^^'^e"^^ + e'^'^ + e"'^'^) . 

However, the function f{d,c,P) does depend on c, as can be easily verified (contradiction). 

The following result shows that, for a given OFBM, one can take the same parameter A in 
the spectral representation (1.3) for all exponents H of the OFBM in question. 

Theorem 6.1 Let Bh be an OFBM having the spectral representation (1.3). If H\,Hjj G £{Bh) 
and Ax, A^j are the two matrix parameters in (1.3) associated with Hx, Hj^, respectively, then 

AxA*x = Ar,A;. (6.7) 

In particular, one may choose the same matrix parameter A in (1-3) for every choice of H & 
£{Bh). 
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Proof: It is enough to show (6.7) with a commuting exponent := Hq (see (5.10)) and the 
associated matrix Ajj := Aq. For simplicity, let H = Hx, A = Ax- We know that 

H - Ho = D - Do=: Ae WCqW'^ = T{Gh), 

where Co C so{n). We can thus write A = WLW~^ with L G so{n). The uniqueness of the 
spectral density of OFBM implies that, for a; > 0, 

x-^AA*x-'^'' = x-^'>AoA*x-^o^ 

or 

^-(Do+A)^^*^-(Do+A)* ^ x-''°AoA*X-^o. 
Since Dq is a commuting exponent and A G T{Gh), then Dq and A commute. Hence, 

x~^°x~^AA*x-^*x'^'o = x-^°AoA*oX-^o 

and x-^AA*x-^'' = AoA*q, i.e., x-^W~^AA*W'^x^ = W'^ AoA*qW-'^ . By differentiating 
with respect to x, we further obtain L{W-^AA*W-^) = (W-^AA*W-^)L, that is, L and 
W-^AA*W-^ commute. Then W-'^AA*W-'^ = W-^AoA*oW-^ or AA* = AqA^. The last 
statement of the theorem follows from (6.7). □ 



A Auxiliary results on matrix commutativity 

We begin by proving Corollary 2.1 and Lemma 2.1. The argument draws upon Theorem 2.1. 

Proof of Corollary 2.1: Without loss of generality, we can assume that ji = 1, . . . , jk = k. 
Since M commutes with A, by Theorem 2.1 we can write M = OKO*, where K = diag(Kii, ii'22) 
and Kii G M{k,M.),K22 G M(n - A;,M). The rest of the claim immediately follows. □ 

Proof of Lemma 2.1: Consider the spectral decomposition A = OTO*, O = (oi, . . . , On) G 
0(n) and T is diagonal. Then, since M and A commute, by Corollary 2.1 and the assumption that 
A e we have M = Odiag(Ai, . . . , \n)0* for some Ai, . . . , A„ e R. In particular, M G S{n, M). 
So, assume by contradiction that Xj^ = . . . = Xj^ (possibly k = 1), where Aj^ 7^ Aj for any other 
eigenvalue Aj of M. Without loss of generality, we can assume that ji = = k. But 

since B commutes with M, again by Corollary 2.1 spauj^joi, . . . , Ok} is a A;-dimensional invariant 
subspace of B, which contradicts the assumption. Therefore, Ai = . . . = A„, as claimed. □ 

The next proposition is used in Theorem 4.1 and Remark 5.5. It shows that, for n > 3, the 
group SO{n) is so rich that only a matrix which is a multiple of the identity can contain it in its 
centralizer. For n = 2, one needs to consider instead the entire orthogonal group. 

Proposition A.l Let T G M(n, M). Then T = XI, X e M., if one of the following assumptions 
holds: 

(i) forn = 2, ifC{T)DO{n); 

(ii) for n>3, ifC{r) D SO{n). 
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Proof: We only prove {u). 

First, let n be odd. We have that F commutes with all rotation matrices of the form 

diag(-l,..., -1,1,-1,...,-!), 

(i.e., the eigenvalue 1 is taken in every "position", all the remaining ones being filled in by — 1; 
there is an even number of entries with —1, so the determinant is 1). By Corollary 2.1, this implies 
that the Euclidean vectors ei, . . . , are all eigenvectors of F. Consequently, F is diagonalizable 
and we can write F = diag(Ai, . . . , A„). Since F G M(n, R), then Ai, . . . , G R. 

Furthermore, more generally we have that F commutes with all orthogonal matrices of the 
form 

Pdiag(-1, . . . , -1, 1, -1, . . . , -1)P*, P G 0(n). 

By the same reasoning as above, this implies that all p G ^"""^ are eigenvectors of F. Consequently, 
every vector v G is an eigenvector of F. 

We claim that Ai = . . . = A^. Assume by contradiction that this does not hold. Without loss 
of generality, we can take Ai 7^ A2- Then F(ei + 62) = AiCi + A2e2, but also r(ei +62) = 7(ei + 62) 
for some 7 G M, since 61+62 is an eigenvector of F and F G M(n,M). Thus, (7 — Ai)6i = (A2 — 7)62. 
Since ei and 62 are linearly independent, 7 = Ai and A2 = 7, which contradicts the assumption 
that Ai 7^ A2. 

Now, let n be even. For notational simplicity, take n = 4. F commutes with 
diag(l, 1, -1, -1) G SO{A). Thus, by Theorem 2.1, F has the form 

/ah \ 
c d 

e / ' 
V 9 h) 

all entries being in M. Likewise, F commutes with diag(l, —1, —1, 1) G 50(4). Thus, by Theorem 
2.1, F has the form 

( i 3\ 
_ m n 

o p 
\k I J 

all entries being in M again. Consequently, F must have the form F = diag{a, d,e,h). In other 
words, 61,..., 64 arc eigenvectors of F (and F is diagonalizable). Now, note that F commutes 
with Pdiag(l, 1, -1, -1)P*, Pdiag(l, -1, -1, 1)P* G S'0(4) for all P G 0(4). Thus, we conclude 
that every p G is an eigenvector of F, and thus that every G is an eigenvector of F. 
Consequently, by the same reasoning as for the last argument for the case of odd n, F = A/ for 
some A G M. 

This argument can be extended to the case of general even n by noting that F commutes with 
all rotation matrices of the form 

diag(-l, . . . , -1, 1, 1, -1, . . . , -1), 

where the "window" associated with the pair (1,1) "slides" all the way down the main diagonal. 
Thus, by again considering all conjugacies P G 0(n), we can see that the same argument applies. 
□ 
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Remark A.l The claim can be extended to the case of F G M(n, C) (which is used in Didier and 
Pipiras (2011), Proposition 5.3). When n is odd, we obtain instead that T = diag(Ai, . . . , A„), 
where Ai, . . . , An G C. Moreover, we conclude that every v G (not C^) is an eigenvector. Thus, 
if Ai 7^ A2, we have that 

7(ei + 62) = r(ei + 62) = Aiei + A2e2 

for some 7 G C, which is a contradiction. The proof for the case of even n can be obtained by 
similar adaptation of the argument for the real case. 

As a consequence of Proposition A.l, we obtain the next lemma, which is used in Theorem 

4.1. 

Lemma A.l // 

WiO{n)W^^ = W20W^^, 

where Wi and W2 are positive definite matrices, and O is a subset (subgroup) of 0{n), then 
O = 0{n). 

Proof: Write QO{n)Q-^ = O, where Q := W^'^Wi. Thus, for any Oi G 0(n), O2 := 
QOiQ-^ G 0{n), i.e., {Q0iQ~^){Q0iQ-^y = I. Thus QOiQ-^ = (Q*)-^OiQ*, or, equiva- 
lently, {Q*Q)Oi = Oi{Q*Q). Thus, by Proposition A.l, there exists A > such that Q*Q = XI. 
Therefore, W2 = WiX~^/^, from which the claim follows. □ 

The next lemma is used in Theorem 5.2. 

Lemma A. 2 Let O G 0(3). Assume O has two (real) linearly independent eigenvectors pi, p2, 
both of which are associated with real eigenvalues. Then O has a third (real) eigenvector p^ such 
that p3 ± span^{pi,p2}. 

Proof: It is clear that, under the assumptions, O has only real eigenvalues. We break up the 
argument into subcases. 

(i) Pi) P2 are eigenvectors associated with the same real eigenvalue (without loss of generality, 
assume the latter is 1): if the third eigenvalue of O is 1, then the proof ends. Otherwise, since 
O G 0(3), eigenvectors associated with different eigenvalues are orthogonal. Thus, there 
exists a vector p^ (associated with —1) which is orthogonal to spa.n^{pi,p2}. Moreover, 
since the eigenvalue in question is real and O G 0(3), one can assume that ps is real. 

(ii) pi, p2 are eigenvectors associated with different real eigenvalues (without loss of generality, 
assume such eigenvalues are 1 and —1, respectively): since O G 0(3), pi _L p2. Moreover, 
since O G 0(3), then the eigcnspacc associated with one of the eigenvalues 1 or —1 has 
dimension 2. Assume without loss of generality that such eigenvalue is 1. Thus, there must 
exist another eigenvector p^ associated with 1. Moreover, since O G 0(3), we can assume 
that ps is real. We obtain that ps -L P2- Also, it is clear that one can choose p3 so that 
P3 _L pi, and thus ps _L spanK{pi,p2}- 

□ 



33 



B On the convergence of eigenvalues and eigenvectors 



The next two lemmas are used in Lemmas 4.1 and 4.2. 

The first one, whose proof is straightforward, estabhshes the sense in which the convergence 
of matrices imphes the convergence of the eigenvalues. 

The second one shows that, under more stringent assumptions, the convergence of matrices 
also implies the convergence of the eigenvectors (in the specific sense described in the claim) . 

Lemma B.l Let {Ak}k^^, Aq G M(n, C), and assume that A}^ — )■ Aq. Then, the eigenvalues of 
ylfc converge to those of Aq, i.e., one can form a sequence {(A^, . . . , A^)}^^^ ^ C" of eigenvalues 
of Ak, k & N, such that 

(A^,...,A^)^(Aj,...,A^), (B.l) 
where the vector on the right-hand side of (B.l) consists of eigenvalues of Aq. 

Lemma B.2 Let {Ak}kQfq, Aq G M(n,C) be such that Aq has pairwise distinct eigenvalues. As- 
sume that Ak ^ Aq. Then, for the sequence of eigenvalues in (B.l), there exists a sequence of 
conjugacies {Pk}kef^ ^ GL{n, C) such that, for large k, 

Afc = Pfcdiag(Ai,...,A^)P,-i 

and Pk ^ P for some P G GL(n, C). Moreover, the columns of the limiting matrix P are, in 
fact, eigenvectors of Aq. 

Proof: By Lemma B.l, the eigenvalues of A^ converge to those of Aq. Consequently, for large 
enough k, the eigenvalues of A^, which we can write as A]., j = 1, . . . , n, are pairwise distinct. For 
all J = 1, . . . , n, choose a basis pI,...,Pq G S^"'~^ of eigenvectors of ^o- 

Now for large enough k and for each j, let Vj^^ := {p G C" : A^p = A|,p} be the one-dimensional 
eigenspace associated with the eigenvalue A'^. Vj^k is well-defined for large k because A^, . . . , A^ 
are pairwise distinct. Now let be one eigenvector in ^'^"-i that is closest to p^, i.e., 

Pfe e argminj,gy.^n52"-i P-Po ■ (B-2) 

The compactness of Vj^k fl ^^^-i implies that a solution must exist. 

Fix any j = 1, . . . ,n. Assume that {pj,} does not converge to Pq, i.e., there exists £o > and 
a subsequence {Pf,,} such that 

" >eQ. (B.3) 



Po 



Since {pj,,} C 5^" ^, one can extract a further subsequence {pj,,,} which is convergent, i.e., 

pI„ ^ p" G S^^-\ (B.4) 

Therefore, Ak^'p^.,, = Xl„pl.„, where Ak//pl.„ Aqp", XI„pI.„ XqP"- Thus, p" G S^^-i is an 
eigenvector of associated with the eigenvalue Aq. Then, since the eigenvalues of are pairwise 
distinct, we can write p" = e'^p;J for some 9 G [—it, it). Note that e~*^p^„ is an eigenvector of A^ 
associated with the eigenvalue A^^.. Thus, from (B.2), (B.3) and (B.4), we obtain 



< £o < 



= e 



Pk" 



0, k" GO 



(contradiction). Therefore, pi 



Pq, j 



1, . . . , n, and we can define Pk = (pj^, ■ 



^Pk) 



Since 



such eigenvectors must be linearly independent, we can write Ak = Pfediag(A^, . . . , A^)P^ ^, and 
Pk^P, where P = {pI...,p1^). □ 
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C Additional results on matrix representations 

The results in this section are used in Section 4. The proof of Lemma C.l is omitted because the 
argument is standard. 

Lemma C.l Let M be a matrix in M{n,M) and let pi,...,pk be linearly independent vec- 
tors in M". Assume span^{pi, . . . ,pk} is an invariant subspace of M. For pk+i, ■ ■ ■ ,Pn 
such that pi, . . . ,pk,Pk+i, ■ ■ ■ ,Pn is a basis of M", define the conjugacy matrix P = 
(pi, . . . ,Pk,Pk+ij ■ ■ ■ iPn) € GL(n, M). Then we obtain the representation 

V M22 y 

where Mil GM(A:,M),M22 G M(n - A:, M), M12 GM(fc,n-fc,M). 

Corollary C.l Let S G S{n, M). // the real orthonormal vectors oi, . . . , generate an invariant 
subspace of S, then 

where O := {oi, . . . ,Ok,Ok+i, ■ ■ ■ ,On) G 0{n), Su G S{k,W),S22 G S{n — k,W). An analogous 
conclusion holds if we replace S G S{n,M.) with L G so{n). 

Corollary C.2 Let S G 5(n,M). Each k-dimensional real invariant subspace of S is generated by 
a set of k real eigenvectors of S. Moreover, the orthogonal subspace is also a {n — k)- dimensional 
invariant subspace. 

Proof: Assume that span]g{oi, . . . , Ok} is a A;-dimensional real invariant subspace of 5, where, 
by the Gram-Schmidt algorithm, we can assume that the real vectors oi , . . . , Ofc are orthonormal. 

Then the representation (C.l) holds for S, where the first k vectors of O are oi, . . . ,0^. Consider 
the spectral decompositions Su = Piidiag(Ai, . . . , Ayk)Pj*]^ and ^22 = -f22diag(Afc_|_i, . . . , A„)P22- 
We can write 

S = Odiag(Pii, P22)diag(Ai, . . . , A^, A^+i, . . . , A„)diag(PA, P2*2)0*- 

Thus, the column vectors in the matrix (oi, . . . , 0k)Pi\ are (linearly independent) eigenvectors of 
S. Moreover, since they are all linear combinations of the vectors oi, . . . , o^, 

spanR{(oi, . . . , Ofc)Pii} C span„{oi, . . . , Ofc}. (C.2) 

On the other hand, since Pn has full rank, then equality holds in (C.2). Thus, the first part of 
the claim is proved. 

The second part can be shown in a similar fashion. □ 
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